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Abstract 

We obtain an effective parametrization of the bulk electronic struc- 
ture of InP within the Tight Binding scheme. Using these parameters, 
we calculate the electronic structure of InP clusters with the size rang- 
ing upto 7.5 nm. The calculated variations in the electronic structure 
as a function of the cluster size is found to be in excellent agree- 
ment with experimental results over the entire range of sizes, estab- 
lishing the effectiveness and transferability of the obtained parameter 
strengths. 
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Effect of quantum confinement on electronic properties of semiconduc- 
tors has been a subject of intense study in the recent years [T] due to their 
potential applications in various devices [2] . In the nanometer regime, where 
the extensions of the electron and hole wave functions are constrained by the 
particle size, the electronic and optical properties are drastically different 
from the bulk material. For example, the band gap of a semiconductor such 
as CdS can be varied from its bulk value of 2.5 eV to about 4.5 eV for ~ 2 nm 
sized clusters ]3] . This size quantization effect can be qualitatively explained 
by the simplest quantum mechanical particle-in-a-box problem where, upon 
decreasing the size of the box, the energy separation between various levels 
increases due to the confinement. Quantitatively, the effective mass approx- 
imation (EMA) [3] describes such a size dependent energy shift. The EM A, 
though suitable for larger clusters, fails to correctly predict the band gap vari- 
ation with the size for smaller clusters, grossly overestimating the band gap of 
the small sized semiconductor nanoparticles. Starting from the atomic limit, 
the ab initio approaches, though accurate and hence desirable, prove to be 
computationally expensive even for small sized ( ~ 1.5 nm) clusters [5]. An 
intermediate approach that uses screened atomic pseudopotentials obtained 
from self-consistent first principle calculations has been applied to II- VI semi- 
conductors H and recently to InP [7J |H] • The results are encouraging but 
this technique too becomes computationally demanding for larger clusters; 
moreover, such approaches are also somewhat empirical in nature in contrast 
to the first principle calculations. Since quantum confinement effects are ap- 
preciable for cluster sizes comparable to the excitonic diameter in the bulk, 
which can be as large as 15 nm in InP, it is obvious that an alternate method- 
ology is required to treat this problem. A semi-empirical tight-binding (TB) 
scheme that borrows parameters from the bulk band structure is less time 
consuming and has been applied with reasonable success to II- VI [9 , III- 
V [TUJ and group IV [TUl HI] semiconductors though the transferability of 
bulk parameters to the nanometer regime is a controversial issue. Our recent 
attempts O EH] in this direction suggest that some of the difficulties arise 
from an improper or inaccurate tight binding (TB) parametrization of the 
bulk band structures. It has been observed that the nearest neighbor only 
model cannot describe the band structure accurately and that inclusion of 
the next nearest neighbor interactions in the basis set is necessary to obtain 
a good description of the bulk band dispersions. [T3] In the present letter, we 
show for InP that a proper parametrization that includes the next nearest 
neighbor interactions leads to an excellent description of experimental results 
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over the entire range of sizes. 

Following our recent successful approach, ^3] we employ a TB method 
with the sp 3 orbital basis on In atom and the sp 3 d 5 orbital basis on the P 
atom. Along with the nearest neighbor In-P interactions, our model also 
includes the second nearest neighbor, In-In and P-P interactions. The s* 
orbital used in the earlier TB models [HI EH El is not required in the present 
approach. 

The TB Hamiltonian is given by 

H = E e u4iaa + E XX*W°Ji' + h.c.) (1) 

il ij W 

where a\ and an are respectively the creation and annihilation operators for 
electrons at the atomic site, % in the I th orbital. The onsite energy for the 
orbital I at the site i is given by e#. The hopping interaction strengths tfj de- 
pend on the type of orbital and geometry of the lattice and are parametrized 
using the Slater-Koster scheme . In order to obtain the best estimates for 
the onsite energies and the Slater-Koster parameters, we obtained the band 
dispersions in InP along various symmetry directions using the Linearized 
Muffin Tin Orbitals within the Atomic Sphere Approximation (LMTO— ASA), 
after converging the calculations with 43 k points in the irreducible part of 
the Brillouin Zone. ^1 The band gap obtained from the LMTO method 
is 0.63 eV; an underestimation of the experimental value of 1.4 eV. Thus, 
we correct the band gap to match the experimental value by shifting the 
unoccupied bands by the appropriate amount. Then, we fit these ab-initio 
band dispersions with those obtained from the TB model within a least- 
squared-error approach by varying the TB parameter values. Fig. ^ shows 
the comparison between LMTO band dispersions and the TB fitted band dis- 
persions with optimized parameter values (Table I), exhibiting an excellent 
agreement over the entire valence and conduction bands. It should be noted 
here that it is not sufficient to get only the energies at various symmetry 
points accurately described, as has been emphasized in the past. The effec- 
tive electron mass from our LMTO and TB calculations, deduced from the 
curvature of the lowest conduction band, are 0.084 and 0.09, respectively, 
showing an excellent fit of the curvature of the band dispersions near the 
T point. They are also in good agreement with the reported experimental 
values ranging between 0.068 to 0.084. ^7] In order to describe the evolution 
of the electronic structure of clusters with the size, it is absolutely essential 
to obtain a reliable description of the curvature of the band dispersion at 
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the extremal points. The reliability of the present model and the parameters 
obtained arises from an accurate description of both the energies and the 
curvatures. [TH] 

Similar to our previous studies of Mn-doped GaAs nanoparticles, [TH] 
the clusters we generate consist of a central In atom surrounded by the four 
nearest neighbor P atoms, followed progressively by shells of In and P atoms, 
respectively. The effective diameter of the cluster is calculated assuming that 
the particles are spherical in shape using the formula 

d = a[*£] ^ (2) 

where a is the bulk lattice parameter (5.861 A for InP [20 ) and N at is the 
number of atoms in the cluster. In our calculations, the largest cluster has 
N at = 9527 atoms and a diameter of 77.1 A containing 60328 orbitals. We 
passivate the clusters with hydrogen atoms at the outermost layer to remove 
the dangling bonds and obtain the eigen-value spectrum for clusters with 
different sizes using the Lanczos algorithm 21. 

Variations of energies of the top of the valence band (TVB) and the bot- 
tom of the conduction band (BCB) as a function of cluster diameter are 
shown in Fig. |21 The results clearly show the systematic shifts of both TVB 
and BCB away from those of the bulk, marked by dashed lines, with decreas- 
ing cluster size as a consequence of quantum confinement. The positions of 
TVB and BCB for four specific cluster sizes have been reported in Ref. [8]; 
we show these results in the same figure indicating a good agreement for 
the common range of sizes. [2H| In order to compare the present results with 
the experimental estimates of band edges obtained from the optical absorp- 
tion experiments, it is necessary to take into account the excitonic binding 
energies at different cluster sizes; we use the form [22J 3.572 *e/ed for the ex- 
citonic binding energy, where d is the diameter of the cluster and e is the the 
dielectric constant, 12.4 for InP. Thus derived calculated band-gap variation 
with cluster size is shown in the inset in Fig. 3, exhibiting a systematic and 
pronounced variation. The dependence of the band gap on the diameter of 
the cluster is expected to be given by an inverse square, l/d 2 , law according 
to EMA. However, such a relationship is completely incompatible with the 
present results. We find that l/d x with x = 1.04 provides the best fit of the 
calculated results, as shown by the dotted line in the inset, suggesting an 
approximately l/d dependence, instead of the l/d 2 one. However, this fit is 
not satisfactory, particularly for large cluster sizes (see inset, Fig- EI) - We find 
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that the expression, 100*(5.8c/ 2 + 27.2<i + 10.4)~ 1 , provides a good description 
of the band gap variation with the cluster diameter over the entire range of 
sizes, as shown by the solid line in the inset. 

The variation of the band gap, AE g , with respect to the bulk band gap 
is compared with experimental values and with results obtained from other 
calculations in the main frame of Fig. El The excellent agreement of the 
present calculated results with the experimental data (221 12H 123 Mj\ over 
the entire range of sizes is evident from this comparison. The expected band 
gap variation on the basis of EMA [27], shown in the figure with a thin 
line, is in gross disagreement with the experimental results. Though the 
parameterized pseudopotential result jS] qualitatively predicts the variation 
of band gap accurately and tends to agree with the experimental band gap as 
one approaches the atomic limit and the bulk limit, it slightly underestimates 
the band gap in the nanometer regime. The sp 3 d 5 s* parametrization scheme 
that was recently reported [TUj overestimates the band gap systematically 
over the entire range of sizes [2*%] . 

In conclusion, we have shown that an effective parametrization within 
the tight binding model allows us to describe the bulk electronic structure of 
InP without the need to invoke the fictitious s* orbital. We have also estab- 
lished the transferability of these parameters to obtain electronic structures 
of clusters over a wide range of sizes; thus calculated variation in the band 
gap with the cluster size is in excellent agreement with experimental results. 
Taking into account the accuracy of this method, the ease of implementation 
and the computational cost, it should be possible to apply this method to 
perform molecular dynamics of large clusters. 

We thank O. K. Aderssen and O. Jepsen for the LMTO codes, P. Mahade- 
van for the LAPW calculations and the Department of Science and Technol- 
ogy for funding the project. 
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The results of Ref. [Hj were rigidly shifted in order to ensure that the 
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order to be consistent with the present convention. 
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Figure 1: Band structure obtained from LMTO and fitted using TB 
parametrization for bulk InP. The zero of the energy corresponds to the 
top of the valence band. 



Figure 2: Variation of top of the valence band (TVB) and bottom of the 
conduction band (BCB) with effective diameter of the clusters. The pseu- 
dopotential results from Ref. [H] are also shown. The dashed lines show the 
TVB and the BCB for bulk InP. 



Figure 3: Comparison of calculated band gap variation with experimen- 
tal data and other theoretical results. The inset shows the 100 * (5.8c/ 2 + 
27.2c/ + 10.4)- 1 fit (solid line) and the l/d im fit (dotted line) to the present 
calculations. 
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Table I 



TB parameters (in eV) obtained from least squared error fit to LMTO band 
dispersions for InP and parameters for hydrogen passivation. 
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